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The production of heavy quarkonium as QQ bound-states in hadron-hadron collisions is consid¬ 
ered within the framework of a space-time description, combining parton-cascade evolution with a 
coalescence model for bound-state formation. The ‘hard’ production of the initial QQ, directly or 
via gluon fragmentation and including both color-singlet and color-octet contributions, is calculated 
from the PQCD cross-sections. The subsequent development of the QQ system is described within 
a space-time generalization of the DGLAP parton-evolution formalism in position- and momentum- 
space. The actual formation of the bound-states is accomplished through overlap of the QQ pair and 
a spectrum of quarkonium wave-functions. This coalescence can only occur after sufficent gluon radi¬ 
ation reduces the QQ relative velocity to a value commensurate with the non-relativistic kinematics 
of these bound systems. The presence of gluon participants in the cascade then is both necessary 
and leads to the natural inclusion of both color-singlet and color-octet mechanisms. The application 
of this approach to pp (pp) collisions from -,/s = 30 GeV - 14 TeV reveals very decent agreement with 
available data from ISR and Tevatron - without the necessity of introducing ht parameters. More¬ 
over, production probabilities are calculated for a complete spectrum of charmonium and bottonium 
states, with the relative significance compared to open charm (bottom) production. An analysis of 
the space-time development is carried through which sheds light on the relevance of gluon radiation 
and color-structure, suggesting a correponding experimental investigation. 


PACS numbers: 14.40.Gx, 13.87.Fh 





I. INTRODUCTION 


The experimental analysis of quarkonia production in high-energy particle collisions has brought in the last couple of 
years exciting discoveries, both vap + p collisions at the Fermilab Tevatron Q , and in Pb + Pb collisions at the CERN 
SPS ||]. In the first case, the interpretation of the Tevatron data of quarkonium production in p-|-p at ^/s = 1800 
GeV, revealed a drastic disagreement with expectations from theoretical calculations that described lower-energy data 
very well. Not only was the quarkonium yield underestimated by a far more than an order of magnitude (in particular 
for J/ip and ip'), but also the observed p_L-dependence of the cross-sections was very different from what was predicted 
earlier. This has become known as the ‘CDF anomaly’. In the second case, a significant suppression oi J/ip production 
observed in the CERN heavy-ion data 0 oi Pb + Pb at Ebeam = 158 A GeV {^/s = 17 A GeV), implied that a naive 
extrapolation from experiments with lighter nuclei substantially overestimated the J/ip production. Clearly, these two 
discoveries concerning quarkonium production, namely, ‘enhancement’ in p + p, and ‘suppression’ in Pb + Pb have 
rather different origin. However, if one is to comprehend the underlying mechanisms in ion-ion applications from 
anything close to fundamental QCD, then both aspects have to be considered. In the case oi p + p (and also e+e“) 
vital theoretical progress based on QCD’s hrst principles has been made in the last few years. In the case of 
heavy-ion collisions the degree of theoretical insight is still perhaps more suspect, at present mainly covering ad-hoc 
prescriptions or somewhat oversimplihed models, however with promising exceptions M. 

To set the stage, I recall that the recent achievements in understanding quarkonium production in hadronic collisions 
at high energy and momentum transfer were essentially stimulated by two ‘observations’ inferred from Tevatron data: 
First, gluon fragmentation (and not direct QQ production) appears to be the dominant source of prompt quarkonia 
at high energies, i.e., a process that involves the high-p_L production of a virtual gluon, followed by its subsequent 
evolution and decay into QQ. Second, a substantial fraction of quarkonium production at very high energy arises 
from intermediate color-octet fluctuations during the evolution into a final-state color-singlet QQ pair, a contribution 
which is small at ISR energies but becomes very prominent at Tevatron energy and beyond. 


Previous to these observations, calculations of charmonium production were limited to the use of the so-called color 
singlet model R-jl] (for a recent review of heavy quarkonium production, see jM). A theoretical breakthrough was 
achieved by Bodwin, Braaten, and Lepage |11|, who developed a rigorous QCD formalism by marrying Perturbative 
QCD (PQCD) and an effective field theory within Non-relativistic QCD (NRQCD). These authors combined in sys¬ 
tematic manner the short-distance, relativistic physics of QQ production and evolution, with the non-perturbative. 
supposedly non-relativistic dynamics of the formation of hadronic QQ bound-states. This novel theory |12[] exploited 
some unique features of the quarkonium bound-states. Firstly, the high-pj_ production of a QQ pair, directly by a 
hard parton collision or indirectly via gluon fragmentation, is a truly perturbative process with the heavy quark mass 
Mq naturally bounding the typical momentum scale from below. Secondly, the internal QQ dynamics is inherently 
non-relativistic, since in order to have a finite probability for producing a QQ bound-state, the pair necessarily must 
have very small relative momentum, despite potentially large total momentum. 


Quantitatively, an important ingredient is the presence in quarkonium physics of several energy (or momentum) 
scales, well separated if the relative velocity v of the QQ satisfies |v| = w <C 1. The three important scales ||^,^ are 
(i) the heavy mass Mq , setting the lower scale for the relativistic QQ production with typical momenta Mq and 
small-size configuratons ^ 1/Mq, (ii) the momentum Mqv, specifying the scale for bound-state formation with typical 
momenta ^ Mqv and mean bound-state size ^ 1/Mqv, and (iii) the kinetic energy Mqv'^, which if small ensures that 
a nonrelativistic treatment of the bound-state may be reasonable mi]- If u <C 1 then Mq ^ Mqv Mqv'^ and 
a double series expansion in the strong coupling as and the velocity parameter v is meaningful. Since u cx 1/ In Mq , 
typical values for the heavy quark velocities are Vcs = 0.23 and Vbi = 0.08. Bodwin, Braaten,and Lepage 
demonstrated that one may then really use PQCD to compute the various color-singlet and color-octet contributions 
of different QQ Fock state components to quarkonium production. The latter, however, are converted into formation 
probabilities for quarkonium only after multiplication by unknown non-perturbative matrix-elements for the QQ 
conversion into physical states. Nonetheless, it triggered several works addressing te application to e+e“ physics 
and pp collisions ||^,|l^. In the context of aforementionened ‘CDF anomaly’, it was shown that by estimating 
the color-singlet contnbutions from quarkonium potential models, and by htting the color-octet contributions to 
experimental data, the addition of the color-octet mechanism can impressively resolve the disagreement of the p + p 
Tevatron data with previous predictions within the color-singlet model only, n 


* It should be noted that the so-called color evaporation model with an early history 0 , was also capable of including octets 
but is too severely constrained by its assumptions Q to have a useful applicability. 
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• color-singlet production 


• color-octet mechanism 





FIG. 1. Feynman graphs that contribute to XQJ production at 0(0®) through color-singlet and color-octet mechanisms. 
The initial QQ is produced in a Pj state, respectively a o’) state, at short distances and hadronizes into a final-state XQJ 
bound-state at large distances. 


Nonetheless, some critical remarks are at hand here. The dictates of experimental data aside, the physical justi¬ 
fication of the color-octet mechanism remains unsatisfactory. To clarify this statement and to motivate the aims of 
this paper, it is useful to draw a somwhat idealistic picture, to be improved on later in this work, of color-singlet and 
color-octet QQ production in hadron-hadron collisions. In Fig. la a typical color-singlet Feynman graph that yields a 
XQJ is depicted. The 0{al) hard scattering of the incoming partons produces at a short-distance scale an outgoing 

color-less QQ pair in a state plus a gluon. The heavy Q and Q fly out from the initial collision point, almost on 

mass-shell, eventually forming, at large distances, a hadronic XQ J bound-state. In Fig. lb, a color-octet QQ pair in 

a state is created by the hard parton scattering, and the pair is assumed to color-neutralize far away from the 
collision point by emitting an extremely long-wavelength gluon which takes away the color, but virtually no energy or 
momentum. The heavy quark pair then may transform again into a XQ J state. 

Whereas the color-singlet production is clearly plausible, the color-octet mechanism seems to contradict basic 
quantum mechanics. The assumption that an ultra-soft gluon takes away the color, but virtually no energy or 
momentum is just wrong! Such a gluon cannot carry away any quantum_numbers within a finite time. The emitted 
gluon must be sufficiently hard if it is to both neutralize the color-octet QQ and separate from the QQ_as an individual 
and incoherent quantum before the QQ pair hadronizes. Since the typical formation time of a QQ bound-state is 
tform ^/Mqv ^ 1 /to, the gluon must have at least an energy ujg = 0(1 GeV) if it is emitted at early time, and 
even more if it is radiated later. Such a gluon is certainly not a long-distance fluctuation with negligible energy. 
Hence, although the intervention of the color-octet mechanism resolves the Tevatron data paradox, I am inclined to 
pose, provocatively, the question: Can the ad-hoc motivation of the color-octet mechanism not be extended towards 
a more physically intuitive scheme, still consistent with the principles of PQCD and NRQCD? In fact, the probability 
of directly producing sufficiently low relative velocity QQ octet or singlet pairs for appreciable coalescence to occur 
is very small. Thus, in a model which attempts to combine both PQCD and NRQCD to generate quarkonium, 
moderately hard gluon radiation plays an integral part and both singlet and octet pairs enter straightforwardly. The 
gluon radiation serves to reduce the difference momentum to a point where coalescence is more probable, and of course 
automatically enables both singlet and octet mechanisms. 

This is the very point where the work of this paper begins. Inspired by the above discussion, I take steps towards_a 
‘microscopic’ space-time QCD description of heavy quarkonium production. The ‘hard’ production of the initial QQ, 
directly or via gluon fragmentation and including both color^singlet and -octet contributions, is calculated from the 
PQCD cross-sections. The subsequent development of the QQ system is described within a space-time generalization 
of the DGLAP parton-evolution formalism Jl^ , pO[ in position- and momentum-space. The actual formation of the 
bound-states is accomplished through overlap of the QQ pair and a spectrum of quarkonium wave-functions. 

The analysis in this paper is confined to pp (pp) collisions in the energy range ^/s = 30 GeV - 14 TeV, spanned by 
ISR, Tevatron and LHC beam energies. Once the space-time structure of quarkonium production in ‘clean’ pp (pp) 
collisions has been better understood, one may step further to address the observed quarkonium suppression in ‘dirty’ 
heavy-ion collisions, where the space-time information of QQ evolution in and interaction with the nuclear medium is 
most essential. 
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In the approach that I advocate here and hope to present in a convincing way in the remainder, it is very natural 
to incorporate both color-octet and color-singlet mechanisms. A produced cc pair can start in either color state, 
although of course they form charmonium only in the singlet. As I shall show, gluon radiation plays so dominant a 
role in the final coalescence into a bound-state that both singlet and octet states must be involved during the pp or 
pp evolution. At the same time it is of great importance to estimate the unknown non-perturbative effects, embodied 
in the conversion probabilities, from a dynamic coalescence picture |^l|| originally constructed for the production of 
deuterons from neutron and proton constituents of nuclear collisions. 

The very explicit cascade space-time structure, for say a pp collision, permits testing for coalescence after repeated 
gluon emission from a heavy quark pair. The qualitative picture that underlies the present phenomenological approach 
is inspired by the above-mentioned combination of PQCD and NRQCD. Accordingly, I also assume that the genealogy 
of a heavy QQ.bound-state hadron involves three, in space-time well-separated components: first, the production 
of an initial QQ pair (directly or via gluon fragmentation), second, its subsequent evolution into a final-state QQ 
pair, and third, the conversion of the QQ pair into a bound-state hadron Hqq. Imagine a heavy QQ pair emerging 
from a hard process of large momenturn transfer scale, with relative probability given by the appropriate hard QCD 
cross-sections. This initial off-shell QQ pair, color-singlet or color-octet, then propagates in space-time according 
to relativistic cascading, and simultanously, but subject to the uncertainty principle, evolves in momentum space 
according to the DGLAP evolution equations that describe the transmutation to the final-state color-singlet QQ in 
terms of real and virtual gluon emission. Finally, the probability for the QQ pair to form a hadronic bound-state 
Hqq is assumed to be given by the overlap of the Q qnd Q wave-functions with the bound-state wave-function. I will 
discuss and examine each of these three, essentially factorizable, building blocks of the model in detail. 

This novel approach to quarkonium production and bound-state formation brings improvements on three distinct 
levels when compared to previous work based on combining (lowest order) PQCD matrix-elements with heavy quark 
effective theory in NRQCD: 

1. A dynamical cascade description is employed that allows to trace the detailed history of the evolution in both space- 

time and momentum space, from the instant of initial QQ production, through the propagation and cascading 
of the pair and its gluonic offspring, and up to the formation of the final bound-state hadron. 

This is to be seen in strong contrast to the previous calculations, which commonly were of a static nature, in 
the sense that no reference was made to the microscopic time-development of the quarkonium system. Nor were 
attempts made to account for ‘real’ gluon emission or possible rescattering of the quark-antiquark composite 
during its evolution. In the present approach on the other hand, the evolution is governed by the dynamics 
itself, including both ‘deterministic’ gluon emission proceses and ‘statistical’ scattering or absorption processes. 
The space-time history is followed causally and provides at any given time the trajectories, coordinates and 
interaction vertices of the particles. 

2. The PQCD aspects of the quarkonia were in earlier work described in terms of lowest order perturbative matrix- 

elements for the production of a quark-antiquark pair in a given state of angular momentum, spin and color-singlet 
or color-octet [^0 ). It was assumed that the initially produced QQ and the final bound-state forming QQ were 
identical states during the perturbative evolution, except for the (in some works) included momentum scale 
dependence. Thus, there was no chance for the initial pair to change its spin or color state after being produced 
by the emission of one or meore perturbative gluons. 

In contrast, the present approach allows the initial QQ to be produced statistically in any spin and color 
combination and then to change its state by gluon emission, absorption or scattering. Thus, the final QQ 
quantum state can be altered and the eventual yield of bound-state hadrons is a statistical outcome of the 
cascade development. Further, matrix-element calculations, though in principle exact, are in practice limited to 
first and second order in PQCD, suffering from the technical complications that enter with increasing order in 
perturbation theory. The cascade description, on the other hand, is based on the DCLAP equations. In principle 
the entire series expansion in Us is included by resummation, although admittedly in an approximate fashion 
and limited to the leading log approximation. Many observables in high-energy e+e“ or pp collisions, for which 
even second order matrix-elements calculations are insufficient, are reproduced with high accuracy in the parton 
cascade description p^ . 

3. The probability for the non-perturbative conversion of the final state QQ state into a bound-state hadron Hqq with 

specific total angular momentum J, orbital angular momentum L and total spin S has been modeled throughout 
the literature by simply lumping all non-perturbative dynamics into the value of the bound-state (radial) wave- 
function TZo{r) for L = 0 (or its derivative 7?.((r) for L = 1) at zero QQ separation r = jrg — rgl =0. The 
underlying motivation for this approximation is borrowed from positronium annihilation, assuming that the 
bound-state hadron emerges from the annihilation of Q and Q, and that the probability for QQ —> Hqq may 
be approximated by the wave-function behavior uniquely at (or near) the origin. Not only does this approach 
neglect the finite width of the bound-state wave-function, but also implies that the probability given by TZq 
or TZ'i is independent of time, hence independent of the actual spatial QQ configuration along time-dependent 
trajectories for the Q and Q. 

In the present approach these shortcomings are absent, since it enables to evaluate at any given time the 
probability for QQ —> Hqq as a function of both relative sparations r and relative momentum k of the Q and Q, 
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using the full expression for the wave-function overlap integral (admittedly, model wave-function), rather than 
simply rely on the value of the wave-function at zero r._ This might be termed a dynamic coalescence picture. 
Moreover, I do not interpret the conversion process as QQ annihilation, but as a coalescence process, analogously 
to, e.g., deuteron production in nuclear physics. Simple but effective techniques have been developed and tested 
in the literature pi| , which I adopt and apply to the present problem. 

The remainder of the paper is organized as follows. Sec. II is devoted to introducing the theoretical framework 
corresponding to the qualitative picture drawn above. Based on the factorization assumption of short-distance PQCD 
production and evolution of heavy QQ pairs, and the long-distance coalescence process with bound-state formation, the 
total squared amplitude for quarkonium production in hadronic collisions is obtained as the product of three functions 
corresponding to QQ production, parton cascade evolution, and coalescence, repectively. It is outlined how the 
resulting quarkonium cross-section can be calulated by simulating the space-time development of the complete particle 
system of a collision event with Monte Carlo methods as implemented in the computer program VNI [^. In Sec. 3 the 
formalism is then specifically applied to pp {pp) collisions in the energy range -^/s = 30 GeV - 14 TeV, corresponding 
to center-of-mass energies between ISR, Tevatron, and LHC experiments. Comparison with experimental data from 
ISR and Tevatron for both charmonium and bottonium production is made, and the underlying dynamics specific 
to the model calculations, is examined in detail, in particular with respect to the relevance of gluon radiation, the 
contributions of color-singlet and color-octet QQ systems, and the characteristics of the space-time structure. Sec. 4 
is reserved for concluding remarks and a brief discussion of future perspectives of the model, especially with respect 
to application to heavy-ion collisions and the issue of quarkonium suppression in dense nuclear matter. 


II. HEAVY QUARK PRODUCTION, EVOLUTION & COALESCENCE TO BOUND-STATES 


Focus of interest is the production of quarkonium states Hqq with definite quantum nubers J,L,S on both the 
charm and bottom sector, for which the following nomenclature is employed: 
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( 1 ) 


Following the outline of Sec. 1, the total amplitude for the production of a quarkonium hadron Hqq in a collision 
of beam particles A, B may be written as a factorized ansatz. As schematically illustrated in Fig. 2, the relevant 
amplitude can then be represented as a product of three components: 

(i) the short-distance perturbative amplitude for production of a time-like fluctuation, being either an off-shell heavy 
pair c* = QQ or a virtual gluon c* = g, plus associated initial-state gluon- and gg-production in a hard collision 
of partons in the incoming beam particles A, B (e.g. p,p) with large momentum transfer Q^ ^ A^; 


(ii) the perturbative evolution of the produced fluctuation c* into a final-state color-singlet [QQY^^ pair with 
definite total spin S, relative and total angular momentum L and J, either by QQ —> [QQY^^^^'> {direct 
production), or by g ^ [QQY’"^{gluon fragmentation), in both cases accompanied by emission of other 
gluons and qq pairs; 

(ill) the long-distance non-perturbative amplitude for the conversion of the final-state color-singlet [QQY^^ as¬ 
sumed to occur by coalescence, into a bound-state hadron HqqY^^^Lj), as categorized by (||). 


Hence, the formula for this factorized ansatz can be expressed as 

A{AB ^ Hqq+ X) = A{AB ^ c +ng + n q + 
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where ng,n'q,n'q etc., are the additional gluons and anti-quarks, respectively, being emitted from the incoming and 
outgoing partons of the hard vertex, and during the cascade evolution (c.f. Fig.2). 


• QQ directly from hard vertex 


A 


B 
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• QQ via gluon fragmentation 


A 


B 




FIG. 2. Schematics _of' quarkonium production in the present approach: Anjnitial QQ pair is produced by a hard collision 
afo —> c*d where c* = QQ (direct production) or c* = g with subsequent g —> QQ (gluon fragmentation). The PQCD evolution 
includes real and virtual radiative processes off the incoming and outgoing partons. The coalescence of the resulting hnal QQ 
pair yields a bound-state hadron with specific quantum numbers J, L, S. 




In order to compute physical observables such as the total production cross-section, one needs to calculate the square 
of the amplitude . On account of the above factorization assumption, interference between the three subamplitudes 
is ignored, and therefore the total squared amplitude becomes a product of three squared subamplitudes, each of which 
may be interpreted probabilistically. Henceforth, I denote these squared subamplitudes by 

(i) 7^^. for the production process resulting in c* = QQ, or c* = g with g —s- QQ, 
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(ii) 5^. for the perturbative evolution in space-time and momentum space, 

(iii) Cjj^g for the coalecence mechanism and formation of a bound-state Lj). 

Indicating in their arguments the relevant momentum and space-time variables, one can therefore express the square 
of as follows: 

\A{AB ^ Hqq+X)'\ = |fci,A:2,fi,r2,t)|^ 

= X! ^c* (M^Q^s|R■o,^o) X I fci,A:2,ri,r2,t) x | fci, fc2, ri,f2, t) . (3) 

c*=g,QQ 

The quantity characterizes the hard scale at the production vertex, = 0(1-2 GeV) is the mass scale that marks 
the boundary between perturbative and non-perturbative regimes, and ^/s is the center-of-mass energy of the colliding 
beam particles A, B. The four-momenta and mean trajectories of the Q and Q, respectively, are defined by 

K = (fc°,k,) rf = (r°,r,) = (ri(t) ) , (4) 

where t = 1, 2 for Q and Q, respectively, and, 

t = \{ri+r 2 ) Ro = (5) 

is the time t, referring to the specific Lorentz frame in which one wishes to describe the dynamical evolution, while to 
and Ro labels the point of time and mean position at which the hard collision occurs. 

Separate discussions of the three functions V^, and are assigned to the next three subsections. The total 
production cross-section for quarkonium mesons of type Hqq is then obtained, in Section IID, by summing over all 
contributing channels and multiplying (^) with the proper flux- and phase-space factors. 


A. Production 


The square of the production amplitudes A {AB ^g* + ng + n'q + (H) gives the probability for 

materializing a time-like fluctuation c* = QQ or c* = g with in a hard collision of partons a,b = at 

the hard scale set by the momentum transfer at the hard vertex. As illustrated in Fig. 3, the squared amplitude 
can be calculated from the usual 2^2 hard-scattering amplitudes for ab c*d, multiplied by the space-like Sudakov 
formfactors Sa.b for the incoming partons a and b and the time-like formfactor Td for the outgoing parton d, resumming 
the contributions of real and virtual gluon radiation. The initial-state partons are weighted with the parton structure 
functions f^/A and fi,/g of the beam particles A,B. Hence, the production probability is 

7^^. (/r^, s I Ro,to) = 


A {AB + + Q)pToducUon 


— ^ ^ J dXa J dX}) f[d ) Sa{Xa')fd | ^o) •) Q I ^o) 

ab 

X |A(a6 ^ c*d)|^ (s,(5^ I Ro,to) x Td(Q^,/i^ | to) • (6) 

where the sum is over all contributing parton channels, and the hard collision ab —> c*d is assumed to occur at time 
to and mean position Ro = ( Ro ). Furthermore, 


s = {PA + PBf g^>lGeV^ = pi . (7) 

are, in that order, the total center-of-mass energy squared, the factorization scale that marks the perturbative boundary 
below which the non-perturbative dynamics is absorbed in the structure-functions, and the momentum transfer scale 
at the hard vertex of the colliding partons a, b, set by the relative transverse momentum of the outgoing parton jets. 
The Mandelstam variables at the parton level, s,),!!, and the lightcone fractions Xa,b of the colliding partons a, b are 
defined as usual 


S = XaXb s t = —XbP±'/se^^ u = —XaP±'/se^^ 


^a,b — 


{kg + kA)a,b 

7^ 


( 8 ) 
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where the incoming parton momenta are parametrized as 


K = (1,0,0,1) k>^ = (1,0,0,-1) , (9) 

while the momenta of the scattered parton c* and of the recoiling parton d, with corresponding jet rapidities j/i and 
2/2 are, respectively, 

kc- = P± (coshyi, l,0,sinhyi) = pj_ (coshy 2 ,-1, 0, sinh 2 / 2 ) . (10) 


• QQ directly from hard vertex 



• QQ via gluon fragnientation 



FIG. 3. Graphical representation of the amplitude for the production and evolution of a QQ pair, analogous to Fig. 2: the 
usual 2^2 hard-scattering amplitudes for ab ^ c*d are multiplied by the space-like Sudakov formfactors Sa,b for the incoming 
partons a and b and the outgoing parton d is dressed by a time-like formfactor Td, which contain the resummmaiMn of real and 
virtual gluon emission. Similarly, the functions Dg, Dq, embody the PQCD evolution of outgoing gluon-QQ system. The 
initial state is weighted with the parton structure-functions fa /a and /b/s of the beam particles A,B. 


The elementary hard process in (D ab —>■ c*d, where a,b = g, q, q, involves the PQCD matrix-elements for producing 
a QQ with definite color and angular momentum quantum numbers to 0(0^). These matrix-elements have been 
calculated for both color-singlet ||,^ and color-octet contributions for a spectrum of QQ states. The explicit, 
rather complicated expressions are well documented in the quoted references. It is however useful to plot in Fig. 4 the 
magnitudes and dependence with momentum transfer of the partonic cross-sections, which are related to the squared 
matrix-elements as usual by da/di = l/(167rs^) \A{ab c*d)\^ The Feynman diagrams that contribute to order 

0(ag) to these partonic cross-sections are depicted in Fig. 5. 
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The factors Xafa/A and Xbh/B in i) give the probablity for hnding the incoming partons a,b in the beam particles 
A, B with momentum fractions Xa,b- The formfactors Sa,b are the probabilities for the incoming partons a, b to emerge 
directly from the structure functions of A, B without additional initial state evolution as illustrated in Fig.3 (conversely, 
1 — Sa,b determine the amount of real and virtual initial-state radiation). The explicit expression for Sa (and similarly 
for Sb) reads 


Sa{xa,tJ.'^,Q'^\t) = exp - 


E 


dt' 


d^r 


dk'^ f^+dz as {{I —z)k'^) 


k'- 


2tt 


la'^a{z) 


f{Xa>,k'‘^) 

f{Xa,k'‘^) 


0 t-t' - 


V 

A:'2 


k' , 


(11) 


where the sum runs over the possible species a' = g, q, q, and the function ^a'^a{z) with energy fraction z = ko a/koc 


is the well-known DGLAP kernel for the energy distribution in the branching o' ^ a. The integration limits 

z~^ = 1 — z~ = 1 — k^/fP of the ^-integration arise from the kinematic constraints of the branching a ^ a'. Last 
but not least, the final-state evolution of the recoiling parton d in Fig. 3 is incorporated in the time-like form factor 


1^ 


Td{Q'^,g? 1 1) = exp 



Us ((1 - 

27r 


{z) 


X 


t-t' - 


V 



( 12 ) 


which is summed over the species d' = g^q^q of the possible daughter species d', and the phase-space boundary for 
the branching is again z~^ = 1 — z~ = 1 — /r^/fc2. 


I emphasize that the expressions dO) and ([l^ are ‘space-time generalizations’ of the standard space-like and time-like 
form factors with so-called angular ordering oerived within the leading-log approximation including coherence effects 
20 2^. By integrating over all times t —> oo and over the spatial variables, one would regain the commonly used 
momentum-space expressions that make no reference to the space-time development of the evolution. Characteristic of 
the above space-time dependent form factors is the appearence of the 6- and d-functions: First, the factor 9{St — ko/k^) 
accounts for the mean formation time tform = ITform = iko/k){l/k) of a gluon emission with the energy ko. By 
enforcing that tform must be smaller than or equal to the considered time interval At = t — t', one ensures that the 
emitted offspring becomes incoherently separate from radiating mother. Notice that instead of d(At — t/orm)? one 
could use a smeared-out distribution such as exp{—At/tform)- The results are rather independent of the specihc form, 
provided only that the choice respects the uncertainty principle. Secondly, the factor (r — vt') with the velocity 
V = k/fco simply reflects the classical propagation of each parton in the cascade evolution, i.e. the spatial development 
is approximated by straight-line trajectories of the partons between branching vertices. 













FIG. 4. 0{al) PQCD cross-sections for the production of QQ states with specihc quantum numbers J, L, S' by a hard parton 

collisions with energy \/i and momentum transfer Q = p±, mediated by gg, qg, and qq collisions. Notice the very different 
magnitude of color-singlet (left) and color-octet cross-sections (right). 


9 













Effective QQ vertices: 
• color-singlet 



• color-octet 



0(q:J) diagrams: color-singlet and color-octet 

• a +9 


=X + X 


5 box graph 
permutations 



I color-octet only ) 


• 9+ g 


Y -g^.grg.'-K 




FIG. 5. Color-singlet and color-octet Feynman diagrams that contribute at 0(0®) to QQ production via the processes (i) 
g + g ^ QQ + g, (ii) g + g ^ QQ + q, and (iii) q + q ^ QQ + g. The top part defines the solid blobs with a fat line attached 
as effective vertices, which are to be understood as insertions in the diagrams of the lower part. 
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B. Evolution 


describes the perturbative evolution subsequent 


The amplitude a(c* [QQ] + m g + m' q + m'q J 

V /evolution _ _ 

to the hard production of the time-like fluctuation c* = QQ or c* = g. If c* = QQ (Fig. 2a, 3a), the QQ pair is 
produced directly at the hard vertex, i.e., without intermediate gluon, so that there is only one stage, namely the 
evolution of the QQ until it becomes a hadron Hqq. If c* = g (Fig. 2b, 3b), there are two stages involved: first, the 
evolution of the time-like gluon through real and virtual gluon radiation until it decays into a (off-shell) QQ pair, and 
second the analogous evolution of the virtual QQ until it converts into a bound-state upon hadronization. Both cases 
can incorporated by representing 5^, as follows: 


evolution 


5^. I fci,A:2,ri,r2,t) = A(c*^[QQY^^’^^^+mg + m'q + mq') 

\ / e 

= Jt J (frc- J ^ J dx J dz {x,z,Q‘^,k^\7c*A) 


zx ^ \zx / 


1 


d: 


Q 


^zx 


(1 — z)x Q \ (1 — z)x 


X2 




(13) 


Here the space-time dependent parton densities D^{x, k'^,fA\ r, t) embody the resummation of the leading-log correc¬ 
tions from virtual and real gluon emsisions to all orders in as according to the DGLAP equations @ , that is, 


^ dt' J d^r j ^ J 


as ((1 - z)k'^) 
- 2 ^^- 


,iz) 


k' , 


(14) 

with ja^ciz) again denoting the DGLAP kernel for branching a ^ c. In ( p^ the function depends on whether 
the QQ pair is produced directly at the hard vertex (fig. 2a, 3a) or via gluon fragmentation (Fig. 2b, 3b): 


{x,z,Q‘^,k^ 


^ ^ f <5(1 - x) 5 (l - (5(r,. - R, 

1 D3{x,Q\k^\r,.,t) 


g^QQ 


iz) 


if c* = QQ 
a c* = g 


(15) 


On the right side, for the case c* = g, the kernel r(z) is now the real branching probability and not the sum of real 
and virtual contributions as the above 7 ( 2 ), since a real QQ pair has to be produced. The two functions are related, 
however, by 


'Ya^biz) = Ta-,bb'iz) - SabS(z - 1) - / dyVa^cc'iv)- 

c,c' ^ 


(16) 


For the case c* = QQ, the product of (5-functions simply cause the integrations over the intermediate c* to collapse 
and enforce the QQ to be produced directly at the mean position R at the hard vertex at time to- 

The physical interpretation of eq. d) in conjunction with ( p^ is rather evident: the product of the functions D 
contain the cumulative effect of the perturbative evolution of the produced c* at Q^ and of the final QQ pair: either (i) 
directly, if c* = QQ, the evolution of the QQ from Q^ to the point when it hadronizes, but at most to the factorization 
point fA in (|^); or (ii) by fragmentation of the gluon, if c* = g, from Q^ to k^, the decay into QQ at k^ and the 
subsequent evolution of the pair from k^ to fP (or above if the QQ form a bound-state beforehand, in which case 
the evolution terminates). 


C. Coalescence 

The final ingredient in the total amplitude (||), is the subamplitude A ( [QQ] for 

the conversion of the final state QQ pair into the bound-state hadron Hqq- Again, I am interested in the square of 
this subamplitude, 
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2 


(17) 


CjLS I fci,fc2,ri,r2,t) 


A 


{[QQ] 


JLS,(1) 




coalescence 


Clearly, in view of the unknown details of the non-perturbative physics, the coalescence of QQ to a bound-state and 
eventually to a physical quarkonium state, requires some model input, and therefore introduces a model dependence 
contrary to the previous elements of the PQCD description of production and evolution. The rationale is here, to 
avoid the introduction of free parameters and develop a transparent and practical bound-state formalism within the 
principles of NRQCD. 


1. Scheme of approximation 


Starting from the amplitude level, the ‘unsquared’ C can be expressed in 
perturbation theory on the basis of the Dyson series 27|, 


the usual fashion of time-dependent 


CjLS I ki,k 2 ,ri,r 2 ,t) 


-i f dt' ) 

J tQ 


(18) 


where IS.E = Eqq — Eh, the wave-functions of the Q and Q are denoted by if, 4’ (tti® brackets [...] indicating the QQ 
coupling to a color-singlet and specific total spin state as described below), and dt is the bound-state wave-function of 
Hqq- The transition operator 

C(Ai^R,^) = C(Ai^R) (19) 


represents the non-perturbative mechanism of conversion from the partonic QQ state to the hadronic bound-state 
Hqq, which in principle may depend on time t, the center-of-mass position R, and of course on the factorization scale 
fj?. Unfortunately the non-peturbative operator V cannot inferred from first principles at present, so that one is forced 
to rely on lattice results, or phenomenological model building. Most previous quarkonium studies simply treated the 
unknown non-perturbative amplitude in (|^) or ( |l8|) as a parameter, i.e., a number being determined either from fits 
to experimental data or taken from NRQCD potential models. 

Instead, I follow a different approach and try to actually calculate approximately this amplitude on the basis of 
the following arguments: It is well known from hadronization studies of jet fragmentation that the mechanism of 
parton-hadron conversion at high energies appears to be a universal and local phenomenon , in the sense that it 

is rather independent of the space-time history of the hadronizing partons and essentially determined by the final-state 
configuration of neighboring partons in a local phase-space cell. Moreover it is known that that the actual details 
of the hadronization dynamics are irrelevant for simple observables as hadron production rates. Hence, it is near at 
hand to assume that the operator V, eq. is approximately of local nature, in the sense that its time variation is 

sufficiently slow that it may be omitted with respect to the QQ motion, and furthermore is at best weakly dependent 
on R (which in the presence of translational invariance is certainly the case). Accordingly, I make the following key 
approximation, pulling V out of the matrix-element (|T^: 


( d- I U(M^R) \[4^) {V) X {^\ [ip^} . (20) 

Returning to (|l^), I insert (|^), take the square, and integrate over the hadronic energy spectrum dEn weighted with 
the level density p{Eh), and obtain so the approximate result for the squared amplitude of eq. (p^: 


C'jLS I ki,k2,Vi,V2,t) 


dEn p{Eh) (V) 


4 

(A^ 



( 4-1 [ip4])f 


2 TTtp{EH) {V) |(VI/| ) 


^ \{'S\[4^)\\ ( 21 ) 

Here I have in the second step exploited the fact that only those final states contribute significantly for which t ~ 
2tt/AE, i.e., AtAE ^ I with At = t — to, and in the last step I have assumed that the energy spectrum p{Eh) oc 1/Eh, 
which is roughly satisfied by the known hadron states. Because the time scale At for the conversion from QQ to Hqq 
is related by the uncertainty principle to the energy difference AE between the quarkonic and hadronic states, the 
approximation (^ ) and (^ip reduces the unknown matrix-element (|^) to the simple overlap matrixelement between 
the heavy quark pair QQ and the bound-state hadron Hqq. 
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8. Wave-functions 


To proceed from (|^), the goal is to model the non-perturbative dynamics of conversion of the final Q and Q 
into the bound-state hadron Hqq by making a specific ansatz for the wave-functions dt and I shall adopt the 

well-established NRQCD bound-state picture for quarkonia 0 in conjuction with coalescence framework developed 
in Refs. |^. As proclaimed before, I restrict the considerations to S-wave and P-wave bound-states with angular 
momentum L = 0 and L = 1, respectively. The corresponding bound-state wave-functions are therefore of scalar, 
respectively vector character in angular momentum space. In the quarkonium bound-state formalism, developed in 
detail in Ref. |Q , the bound-state wave-function is constructed in terms of the overlap amplitude between quark and 
antiquark spinors tp and tp in a given total spin configuration (singlet S' = 0 or triplet S = 1) and the bound-state 
wave-function d' in a given angular momentum state L = 0 or L = 1, to form s state of total angular momentum 
J = 0,1, 2. As usual, I assume that in the restframe of the bound-state the relative momentum |k| = §|k;i — k 2 | of the 
Q and Q with 3-momenta ki and k 2 , respectively, is small as compared to the quark mass Mq. As a consequence, 
one can decompose the (unsquared) amplitude Cjls, eq. (HI). in coordinate space as follows: 


CjLS I fci, fc2,ri,r2,t) 


2Tr6 





-== ^ {LM,SS,\JJ.) vl/^*(ri,r2) 


[tPPJ\g^^\ri,r2) . 


( 22 ) 


Here ( L M, S Sz \ J Jz ) the Clebsch-Gordan coefficient for the coupling of spin S and relative angular momentum L 
to a total J with associated projections S^, M,and Jz, respectively, is the complex conjugate of the bound-state 
wave-function for a L-state with projection M, and 


(l’l,^i;i’2,fc2) = 




Nc 

E 

Cl.C2 = l 


3z, 3j|10)^ ( xS 2|S'5'^ ) V'(l*l,Sl,Ci) V'(i*2 ,S2,C2) (23) 


couples the Q and Q spinors to a color-singlet with total spin S and projection Sz- 

To proceed, one needs to specify the wave-functions 'k and ip, ip. Most plausible would be to use quarkonium 
wave functions from potential model calculations, which however are not available analytic form. Therefore, I choose 
instead to work with harmonic oscillator wave-functions for which it is possible to obtain the overlap integrals involving 
"^,ip,ip, as transparent, closed expressions. Following the coalescence aproach of Refs. [Q, I shall assume that quark 
Q and antiquark Q may be described by wave-packets of width (7^ oc 1/ {kf = — kf), localized in space around 

Ti and in momentum space around k^, multiplied by the usual SU{2) and SU(3) spinors y, ^ for spin and color, 
respectively: 


ip{ri,Si,Ci) = - -exp f- exp(ik,-r*) xl i = 1,2 {= Q, Q) (24) 

Here and in the following the notation of 3-vectors is such that, e.g., refers to the mean position as a classical 
variable, whereas describes the quantum-mechanical fluctuation of a particle’s position around the mean value. 
Notice that the plane-wave factor exp(iki • r^) ensures the uncertainty principle in the sense that a classical particle 
interpretation is only possible if |r| << l/|ki|. 

The bound-state wave function on the other hand, I write as a product of its center-of-mass motion and its 
internal motion, again in terms of Gaussian wave-packets, the size of which generally depends on the particular orbital 
angular momentum state, and of course on the flavor of the heavy-quark pair. Defining center-of-mass and relative 
coordinates and momenta as 




ri + i ’2 


K = ki + k2 


r = ri - r 2 

, ki - k 2 
k = - 


and similarly for the mean positions R, r and momenta K, k, I make the factorized ansatz 

^M(i’i.r 2 ) = 4>^(R) (piiir) . 

The center-of-mass wave-function depends on L in its spatial width, 

I P (R-R)2' 


$^(R) = 


(7rEi)3/4 


exp - 


2 Si 


exp(i K • R) 


(25) 


(26) 


(27) 
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Ylm(j) , 


(28) 


and the relative motion depends on both L and its projection M, 




Ml 


(7rai)3/4 


M exp 


2al 


with the L-dependent prefactors A/) = 1 and Mi = -\/2/3 and the spherical harmonics Y^m depending on the angular 
orientation r = r/r, where r = |r|. 


S. Coalescence probabilities 


Inserting the wave-functions and dt into (|^), integrating over ri = R -|- and ri = R — ^r, summing over 
the angular momentum projections, and then squaring the resulting amplitude Cj^g, one obtains the coalescence 
probability Cjj^g for the conversion of QQ into bound-states Specifically for the S- and P-wave states, one 

obtains the following closed expressions: 

• S'-states (L = 0, J = S' = 0,1): 




JOS 


t] = 




exp 


2 


M5 


(29) 


P-states (L = 1, J = P -|- S = 0,1, 2): 




JIS 




4 1/1 


V2 


Ml 


0^1 —2 
— ^ 

Ml 


—2 


Vi k 


exp 


-dt 
2 * 


Ml 


(30) 


In the above formulae 


and 


I used fc = |k| = 4|ki — k 2 |, r = |r| = |ri — r 2 |, and 


2 _ o 2 I 2 

Mi = 2 cr + aj^ 


Ml 


(31) 


Furthermore, M = ^{cti + cr|), and = M/2, assuming that the 2-body center-of-mass motion of Q and Q is 
unaltered during the coalescence process. 

It is worth noting that the final expressions for the coalescence probabilities in (^^ and (^) are functions of 
the relative QQ momentum and separation only, and do not depend neither on the individual variables k^ and r^, nor 
on the total momentum K or the cm-position R of the pair. This makes the coalescence process a local mechanism, in 
accordance with the aforementioned locality of hadronization in high-energy particle collisions. Also notice that the 
coalescence probabilities are time-dependent, although not explicit in ( p^ and (|^), due to the dynamically changing 
Q and Q momenta and trajectories, i.e., k = k(t) and r = f{t). Finally, the dependence on indicated by the first 
argument of Cj^g, is again implicit in the Q and Q momenta, since per definition p"^ serves as a cut-off for soft gluon 
radiation with energies below p from Q oi Q and hence limits the radiative energy loss of the pair. This in turn 
influences the coalescence to a QQ bound-state, although this dependence is very weak. 

The independent parameters involved in ( p^ and (30) are M (a property of the Q and Q upon coalescence, but 
independent of the bound-states Hqq), and a\ (a property of the bound-states Hqq, depending on the angular 
momentum quantum number L as well as on the ‘flavor’, cc or bb). However, these parameters M and are not for 
free choice. Rather than that, the possible values of M, are given by the invariant (off-shell) masses kf = — kf, 

and hence are dictated by the dynamics of Q and Q itself, whereas aj^ is determined by the Bohr radii of the various 
quarkonium states Specifically, for M = -|- tT|), the widths ai of the Q and Q wave-packets 

are obtained from relating their (time-evolving) invariant masses kf by means of the uncertainty principle to the 
corresponding mean spatial sizes [], 


^ The factors 2/3 in ds^ ) and (p^) arise from the definition of the mean square radii of the heavy quarks and the bound-state, 
respectively, which involve the ratios (r'^) = J dfr'r f{r')/ J dfr' f{Y), where f{r) stands for the Gaussian wave-packets of 
Q, Q and of the bound-state. 


14 









a? = 


3 k. 


2 ’ 


kf = kl, - kj 


> Ml 


(32) 


whereas the widths of the bound-state wave-functions are related to the sizes of the corresponding quarkonium 
hadrons through the measurable electron scattering radii, 

^ ^UHqq) = I^Iq^ Rqq ^ ■ ( 33 ) 

The relevant values of Rqq for charmonium and bottonium states that are used later, are listed in Table 1. 





R -(25+1 

Lj) ( fm ) 




'5o 


'Pi 

"Po 

"Pi 

"P2 

cc 

Vc 

J/4, 


XcO 

XcV 



0.25 

0.27 

0.55 

0.60 

0.63 

0.66 

hi 

rib 

T 

h 

XbQ 

Xbi 

Xh2 


0.10 

0.12 

0.25 

0.29 

0.31 

0.33 


TABLE I. Quarkonium RMS values Rqq used for the widths of bound-state wave functions with dif¬ 
ferent J, L, S in conjunction with the relation (P3|). 


The coalescence probabilities (|2^) and consist (apart from the normalization factors that are always < 1) of two 

Gaussians, one in momentum space, exp(—/2), and one in position space, exp(—r^//i^). The first one ensures the 
conservation of momentum, whereas the second restricts the distance between the coalescing Q and Q to a width which 
is given by the bound-state size, augmented by the uncertainty in the positions of the Q and Q. The essential difference 




between the S-wave coalescence probabilities (|^ and the P-wave counterpart (30), is the factor (a^T^//r‘* -|- v‘ 
present in the latter. For k —s- 0 the P-wave probability Cj^g therefore vanishes, reaches a maximum as r increases, 
and then tends to zero again with an exponential tail. Contrary to this, the S-wave probability Cjgg, for k —> 0, has 
a maximum at r = 0 and damps out exponentially as r increases. This characteristic behavior is plotted in Fig. 6, 
which shows the coalesence probabilities (EM and (BM for the different angular momentum eigenstates as a 


function of r = |r| for various values of A: = |k|. 


4- Relation between Cjls and commonly used quarkonium wave-funetions TZl{0) 


At this point, it is instructive to relate the coalescence probabilities in ( p^ and (|^ with their analogues in 
earlier quarkonium studies, namely the values of the squared radial quarkonium wave-functions at the origin, |Pi(0)p, 
where 


7^L(0) 


drL 


Rdr) 


r—0 


(34) 


with r = |r|. Interrelating |7?.i(0)p and is somewhat ambigous, because the latter coalescence probabilities are 
to be weighted with and summed over all possible QQ configurations with distinct relative momentum and spatial 
separation, while the former |Pi(0)p is just a number, independent of space-time and the momentum distribution of 
Q and Q. Hence, one expects the summed-over coalescence probabilities generally be larger. Nevertheless, one may 
establish as a limiting estimate an approximate correspondance between the Cj^g and the |Pl(0 )P by concentrating 
on the phase-space region where coalescence receives its dominant contribution from, namely at the maximum value. 
This is in accord with the commonly used approximation of using the value of is approximately given by 7?.L(r) at 
r = 0 where the radial wave-function is maximal. 
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Coalescence Probabilities 





r (fin) 




FIG. 6. Typical shape of the coalesence probabilities (j2^ and for the different angnlar momentum eigenstates 
as a function of r = |r| for various values of k = |k|. The values of fiL and ul entering the expressions and are 
determined from an (|33|), using the RMS values of Table 1. 


Consequently, one may write 

where the maximum values of the coalescence probabilities (|^) and (^0[) are 




AL) 


0 for L = 0 
aL for L = 1 ’ 


(35) 


(36) 


and t = tc is the corresponding time argument that can be set to zero, since it is s dummy variable for the present 
execise. The integration is over the typical bound-state size with radius Rqq = Rj)- f^e limit k ^ 0 

one may therefore write for the cases L = 0 and L = 1, respectively. 
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C^j,s{pi^\0,0,0)) |7^o(0)|" R 


2 o3_ 

QQ ’ 


^jis I 0, cn, 0)) 


(fe^O) 




d3_ 

2 -^QQ 


for L = 2 . 


(37) 


Typical values for |72 ,l( 0)P from QCD potential model calculations 


|77i(0)r 

, e.g. for Buctirniiller-type wave-functions, are 


l^o( 0 )|" = 


0.81Get^3 

6.48 Gey3 for 66 


i^iwr 


0.041 GeT43 for cc 
0.057 GeV^ for 66 


(38) 


for S'-wave and P-wave states, respectively, with Me = 1.35 GeV and Mf, =5 GeV. Gomnariiig these values with the 
coalescence probabilities in the limit (|^), one obtains, on account of the explicit expressions (|^) and (|30|), and using 
the values of Rqq(^^^^Lj) from Table 1, the following estimates for S'-states and P-states, respectively: 


1 ^J 05 I 0)0; 0)) ^ J 1.23 for cc Mq C'j^g [y? |0, ai,0)) ^ J 1.15 for cc . . 

Rqq |7^o(0)|" 1.29 for 66 ’ |Pi(0)|" “ \ 1.30 for 66 

One sees that these ratios for ^-wave and P-wave states, respectively, indicate that the coalescence model used in 
this paper gives around 15-30 % larger values for the formation probabilities of quarkonium than the widely used 
approximation based on the values of the radial wave-functions at zero QQ separation. It has to be kept in mind, 
though, that the actual magnitude of the coalescence probabilities can be significantly larger than in ( ^^ when summed 
over all possible QQ momentum configurations and trajectories. 


D. Cross-sections 


In order to make contact with observable quantities such as multiplicities, relative abundances or spectra of produced 
quarkonium mesons, one needs to integrate the space-time dependent total squred amplitude (^ over the complete 
event history of a collision A + B, as well as over all unobserved particle states X. For instance, the total production 
cross-section for quarkonia Lj) is obtained after dividing by th incoming parton flux and setting without 

loss of generality fo = 0 , Rq = 0 , 

da {AB Hqq +X) = da {AB c* X) (s, Q^) x /i^) , (40) 

c*=g,QQ 

where 

da [AB C* X) (s, Q^) = Y fa/A(Xa, M^) fb/B[xb, /i^) Saixa.fM^, | 0) SbiXb, | 0) Td{Q'^, I 0) X 

abd 


2s 


\A[ab c*d)f [s, Q^) ) [2Tr)'^6[ka + h - K* - kd) dxa dxb d'^kc^d'^kd ■ (41) 


and 


(Q^, = J dt J d^rid^r2 J d'^kid^k2 £e-^ {Q^, \ ki, k2,ri,r2,t) (a^^ | ^i, fe, ri, r2, t) . (42) 

Using the kinematic variable definitions (0) -(p^, one can express the triple-differential hadronic cross-section for Hqq 
production in terms of the rapidities yi of the jet c* that initiates the QQ evolution, the rapidity y 2 of the recoiling 
jet d which evolves and fragments by itself, and the common relative transverse momentum Q = p± posessed by each 
of the two outgoing jets c* and d, 

(AB^H^^ + X)= j: ^ d X) . . (43) 


dyidy2dp]_ 

where is the same as above, and 
da 


c*=g,QQ 


dyidy2dp]_ 


[AB^c*X) = ^ a:a/a/A(xa,Ai^) a:&/b/B(a:b,M^) ( ^6^ ^ c*d)\‘^ [s,p'i)] x 

abd ^ ' 

X S'a(xo,A‘^,pi I 0) Sbixb,p^,p'i I 0) Td[p']_,pL^ I 0) . (44) 
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E. Calculation scheme 


Due to its probabilistic character, the above formalism is suitable for implementation in a Monte Carlo simulation 
using the computer program VNI p^, which allows one to trace the microscopic history of the dynamically-evolving 
particle system in space-time and momentum space, from the instant of collision of the beam particles A and B to 
the eventual emergence of quarkonium states plus all other final-state particles. Specifically, following the space- 
time development of the complete system of particles in a collision event, embodies the dynamics of the QQ system 
within the rest of the event, consisting of all other produced partons and the hadronization products, as well as the 
fragmentation of the spectator partons of the beam remants of the original A and B. This allows to account not only 
for the correlations of particles in space, time, color and flavor, but also for the statistical occurence of quarkonium 
formation by QQ coalescence versus open charm/bottom production by recombination of Q and Q with partons of 
the rest of the event. 

Referring for technical details to Ref. I confine myself here to a brief outline of the simulation scheme. The 
time development in position and momentum space of the space-time dependent admixture of partons and emerging 
hadronic excitations, including Q and Q and its final hadronic products, is obtained by following each individual 
particle through its history with the various probabilities and time scales of interactions sampled stochastically from the 
relevant cross-sections. The detailed history of the system can thus be traced by evolving the phase-space distributions 
of particles are evolved in small time steps (At ~ 10“^ fm) and 7-dimensional phase-space d^rd^kdE, from the instant 
of collision between the beam particles A and B, all the way until stable final-state hadrons and other particles are 
left as freely-streaming particles. At any time during the evolution, the state of the QQ system can be extracted, and 
so the final quarkonium yield at the end of of a collisions event. For the case of pp and pp collisions on which I shall 
concentrate in the following, the essential ingredients in this Monte-Carlo procedure are summarized as follows 

(i) The initial state is constructed by decomposing the incoming beam particles, p and p, into their parton substructure 

according to measured structure functions ^ with a spatial distribution given by the Fourier transform of the 
proton elastic form-factor. The so-initialized phase-space densities of (off-shell) partons are then boosted with 
the proper Lorentz factor to the center-of-mass frame of the colliding beam particles. 

(ii) The initial QQ production is obtained by triggering on those hard parton collisions ab c*d that contain a 

c* = QQ or c* = g* with g* —>■ QQ. These hard processes are sampled from the well-documented cross- 
sections , including the direct process and gluon fragmentation process with weights given by the relative 

magnitudes oFthe corresponding cross-sections. 

(iii) The parton cascade development proceeds then by propagating all partons along classical trajectories until they 
interact (radiate, coalesce, or hadronize), including the QQ system as well as te additionally produced partons 
from initial-state space-like radiation off a, b and final-state time-like radiation off the QQ system and the recoiling 
jet d. The production of each parton through the hard scattering or radiation is subject to an individual, 
kinematic-specific formation time Atp^c = 'y/Mp where 1/Mp = l/V^ is the proper decay time of an off-shell 
parton with invariant mass Mp and 7 = E/Mp is the Lorentz factor. 

(iv) The final QQ quarkonium formation is traced by evaluating in each infinitesimal step of the collision event the 
spectrum of coalescence probabilities (|^ and (|3^) from the Q and QQ positions and momenta, and testing the 
values against a random number. At each point m time the total probability determines whether a quarkonium 
state may be produced, and if so, the relative probabilities of the different states give the weights for which 
quarkonium state out of the spectrum is formed. If no coalescence occurs, the Q and Q continue in their 
evolution as indiviual particles as in (iii). 

(v) The emergence of hadrons other than quarkonium involves the recombination of produced partons into pre-hadronic 

cluster, using the Ellis-Geiger scheme in which at any point during the cascade dvelopment any pair of 
nearest-neigbor partons is tested whether their mutual separation favors a cluster formation, or whether they 
proceed in their parton evolution as in (iii). These pre-hadronic clusters decay subsequently into primary hadrons, 
which in turn if excited or resonant states, can decay into final stable particles according to the particle data 
tables |3^ . Again, each newly produced pre-hadronic cluster as well as emerging hadron becomes a ‘real’ particle 
only after a characteristic formation time AQ.h = 'j/Mc^h depending on the invariant mass Mc,h and the energy 
through 7 = E/Mc,h. Before that time has passed, a cluster or hadron is considered as a still virtual object that 
cannot interact incoherently until it has formed according to the uncertainty principle. 


^ The GRV structure function parametrization is used throughout, which describes quite accurately the HERA data even 
at low and very small x. 
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(vi) The beam remnants, being the unscathed remainders of the initial beam particles (in the present case p and p), 
emerge from reassembling all those remnant primary partons that have been spectators without interactions 
throughout the evolution. The recollection of those yields two corresponding beam clusters with definite charge, 
baryon number, energy-momentum and position, as given by the sum of their constituents. These beam clusters 
decay into final-state hadrons which recede along the beam direction at large rapidities of the beam/target 
fragmentation regions. Again, individual formation times of the produced hadrons are accounted for. 

It is important to realize that the spatial density and the momentum distribution of the particles are intimately 
connected: The momentum distribution continuously changes through the interactions and determines how the 
quanta propagate in coordinate space. In turn, the probability for subsequent interactions depends on the resulting 
local particle density. Consequently, the development of the phase-space densities is a complex interplay, which - at a 
given point of time - contains implicitly the complete preceding history of the system. 


Parameter 

Value 

Description 

At 

0.001 fm 

time step by which the evolution of the particle system 

is followed in the hadronic center-of-mass frame. 

j'^QCD 

0.29 GeV 

QCD scale used in and throughout the perturbative 

evolution of the parton distributions. 

1^ 

1.0 GeV 

mass scale below which the perturbative parton evolution 
is cut off, if not terminated before by coalsecence or 

hadron formation. 

M, 

1.35 GeV 

value used for charm quark mass 

Mk 

5.0 GeV 

value used for bottom quark mass 


TABLE II. Parameter values used in the model calculations. 


^ For the case of pp or pp collisions considered in the following, the general variety of particle interactions reduces to only the 
radiative parton poduction associated with the evolution before and after the hard scattering, and the parton-hadron conversion 
through parto-parton recombination or coalescence. 
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III. APPLICATION TO PP{PP) COLLISIONS AT 30 GEV - 14 TEV 


For the production of the charmonium and bottonium S-wave states J/ip and T, respectively, there exists a body 
of experimental data, both from pp collisions at the ISR {■\/s = 30-63 GeV) and from pp experiments at the Tevatron 
(•\/s = 1.8 TeV). Also, the charmonium P-wave states Xc have been observed, although with larger experimental 
uncertainty, since these states are only indirectly measured through their radiative decay into J/ip (unlike the S-wave 
states which due to their narrow resonance character can be observed directly through invariant mass analysis). On 
the other hand, for radial S-wave excitations, such as ip', T', T", data are more scarce, and detailed measurements of 
Xb as well as of the and ^Pi states, rjc, r]b, he, ht, are not reported at all. 

Hence, I will focus primarily on the production of the charmonium states J/ip (^S'l) and Xcj i^Pj), as well as the 
lowest-lying bottonium state T (^S'l). However, since in the model calculations the production patterns of the various 
quarkonium states are intimately entwined by the dynamical evolution, I also discuss somewhat the relevance of the 
other S-wave and P-wave states. 


A. Production cross-sections: postdictions and predictions 


The principal result of this paper is contained in Fig. 7, which collects the quarkonium production cross-sections for 
J/ip, Xc and T, and which is posted here forefront as basis for the subsequent analysis of the underlying dynamics in 
the present approach. Shown in Fig. 7 is the p^-dependence of the cross-sections around central rapidity at 63 GeV, 
1.8 TeV and 14 TeV, as given by eq. (^), 


da 

dp± 


da 

= 2pj_ 2 
y=0 dP± 


for pp{pp) 


yi-i-y2<i 


J/ip + X 

Xc J/'tp + X 
T + X 


(45) 


multiplied by the leptonic branching ratios BjJ/fP = 5.97%, and B{T ) = 2.48%. Confronted 

are the calculated cross-sections (from eqs. (0) and (^^) with the measured spectra from Refs. |^,Q. In order 
to allow direct comparison with experimentalaata, the following quarkonium yields are plotted: The top panel of 
Fig. 7 displays the production of prompt J/ip's only, with the radiative feed-down from as well as from bottom 
decays subtracted. The middle part reflects the production of Xc summed over all J = 0,1,2 states, and implicitely 
obtained through the amount of J /ip's that result exclusively from the radiative decays of Xc with the branching ratios 
Bixco —> J/ip-\-"f) = 6.6-10“^%, B{xci J/V' + t) = 27.3%, B{xc2 —> T/V' + t) = 13.5%. Finally, the bottom 
panel shows the production of T’s, including the radiative feed-down from the Xb states. 

The features of the results in Fig. 7 may be itemized as follows: 


(i) The agreement of the calculated cross-sections for J/ip, Xc and T production with the experimental data at ISR 
and Tevatron is very good. Not only the pj_-shapes of the cross-sections, but also the overall magnitudes are 
postdicted correctly for ISR and Tevatron. The LHG curves may therefore be seen as quantitative predictions. In 
view of the fact that the dynamics of quarkonium production in the model calculations is essentially parameter- 
free (aside from the fixed quarkonium radii in Table 1), this success suggests that the combination of PQCD 
parton evolution with bound-state formation via coalescence may in fact be a reasonable visiualization of the 
underlying physics. 


(ii) The overall pattern of charmonium and bottonium production is quite similar, although with important differences 

in detail. The comparison of the J/ip and T production reveals a ‘flattening’ of the p_L-spectra due to the larger 
quark mass in the latter case (c.f. Table 2), and, connected with this, a much smaller size of the T as compared 
to the J/ip (c.f. Table 1). The confrontation of the Xc and J/ip spectra implies a softer, i.e. ‘steeper’, p± tail in 
the Xc production, which again is plausible due to the larger size of the Xc states. The ‘flattening’ indicates a 

power-behavior in the transverse energy 4 Mq -|- p\ rather than in the transverse momentum p_L itself. 

(iii) Most prominent in the charmonium sector is Xc production, leading also to a strong Xc component in the inclusive 

J/ip yield through radiative feed-down Xc J /'P’ + 7- In fact, the calculations show that more than 30 % of 

the total J/ip yield originates from Xc decays. On the bottonium sector this pattern is repeated, as the model 
calculations indicate. Again, the x&’s provide the largest yield, however, they do not contribute as strongly to 
the T population as in the corresponding case of Xc and J/p. This is due to to the larger quark mass which 
leads to a kinematic suppression at large p±_. 
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B(JAlf dc/dpj(pp /A|/ + X) 



Pj(GeV) 

B(JA\f j do/dp^pp ^ X 



Pj(GeV) 



Pj(GeV) 


FIG. 7. Comparison of calculated and measured cross-sections for quarkonium production around central rapidity at 63 
GeV, 1.8 TeV and 14 TeV. Top part includes prompt J/'tp’s only, with the radiative feed-down from Xc’s as well as from bottom 
decays subtracted. Middle part shows Xc production summed over all J = 0,1,2 states, identified from J/'tp’s that result 
exclusively from the Xc decays. Bottom part contains T’s, including the radiative feed-down from the Xb states. For better 
distinction, the cross-sections have been scaled by factors of 0.1 (63 GeV), 1 (1.8 TeV) and 10 (14 TeV). The full curves are the 
model results. The data points are from Ref. [34] for ISR at 63 GeV and from Ref. [1] for Tevatron at 1.8 TeV. 
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In view of the non-trivial physics of quarkonium production as discussed in Sec. I, the good agreement of the model 
calculations with the available data immediately poses the the question, 

Why does this space-time description, which combines PQCD parton production and evolution with a 
simple coalescence model, agree so well with experiment? 

More specifically, there are three main issues to be examined in the following, namely, 

• Initial QQ production: What is the underlying production pattern, in particular, the relative importance of 
color-singlet and color-octet contributions? 

• Parton evolution: How does the PQCD evolution and the associated gluon radiation affect the quarkonium yield? 

• Coalescence: How does the space-time development of the initial QQ proceed and what are the characteristics 
of the final QQ upon coalescence? 


PP (PP) CLt GeV -14 TeV 



Q (GeV) 



10 

Ecn. (GeV) 


FIG. 8. Top panel: Distribution of values for the pp (pp) energy range ^/s = 30 GeV - 14 TeV, where Q = p± represents 
the momentum juansfer at the hard vertex ab ^ c* d of Fig 2. Bottom panel: Integrated partonic cross-section for the production 
of the initial QQ from the hard vertex ab c* d (directly, or via gluon fragmentation) in th energy range = 30 GeV - 14 
TeV. 
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In the following Subsections these questions are addressed and discussed in detail. Beforehand, I would like to draw 
attention to Fig. 8, which may serve to further some intuition for the beam-energy (-^/s) dependence of quarkonium 
production and the relative importance of contributing partonic production processes. The top part of Fig. 8 displays 
the distribution of = p]_ values for the pp (pp) beam energies 30 GeV - 14 TeV. It reflects the typical momentum 
transfers at the hard vertex ab —> c*d, where a,b,d = g, q, q and c* = QQ (Fig. 2a) or c* = g with g —s- QQ (Fig. 2b). 
As expected, there is a substantial increase in the mean Q^, accompanied by an increasingly prominent power-law 
tail, as one moves from ISR energy 30 GeV up to LHC energy 14 TeV. The immediate consequences are, firstly, an 
enlarging of the kinematic region for producing heavier quarkonium states, and secondly, the increasing importance 
of gluon radiation, the role of which I shall discuss in detail later. The bottom part of Fig. 8 shows the beam-energy 
dependence of the integrated partonic cross-section for the production of the initial QQ from the hard vertex ab ^ c*d 
(directly, or via gluon fragmentation). Three features are self-evident here. First, the cross-section rises by almost two 
orders of magnitude in the range 30 GeV - 14 TeV. Second, the channel ab = qq of light qq annihilation is down by 2 
orders of magnitude, and is therefore completly irrelevant in the energy range considered here. Third, the contributions 
of ab = qg and ab = gg are about equal in magnitude throughout, with the qg-channel slightly dominating. This is 
a very important property that is often overlooked, because one would naively expect that the gg-chtamel dominates 
by far due to the strong presence of low-x gluons in the protons. However, collsions of low-x gluons are suppressed in 
the QQ production, because of the failure to overcome the mass threshold. 


B. Properties of color-singlet versus color-octet contributons 


In addressing the first of the above issues, namely the relative importance of color-singlet versus color-octet contri¬ 
butions to the final-state quarkonium yield, recall from the Sec. I, that historically the ‘color-singlet model’ did 
very well at ISR energies and below, but underpredicts badly the quarkonium productions at the Tevatron oy more 
than an order of magnitude, and thereby suggesting a strong ‘color-octet mechanism’ |^,|l^,|^ contributing at very 
high energies. 

In the present approach, both color-singlet and color-octet production of the initial QQ are self-consistently included, 
by determining the relative probabilities for singlet and octet production from the respective weights of the hard QQ- 
production cross-sections. It is important to stress again that an initial singlet (octet) QQ pair produced by the 
hard parton collision does not necessarily result in a final QQ with the same quantum numbers upon coalescence, 
because the emission of gluons during the evolution changes the color state (and in principle also spin and angular 
momentum) of the pair. In fact, since in the present approach, the final QQ has to be in a color-singlet state in order 
to coalesce to a quarkonium bound-state, this is a crucial necessity to be able to incorporate both singlet and octet 
production of initial QQ pairs. Figs. 9-11 separate the color-singlet and color-octet contributions and discriminate 
their characteristic features. 

Fig. 9 compares the relative proportions of the various included charmonium states at Tevatron energy ^/s = 1.8 
TeV (the pattern is very similar on the bottonium sector, not shown here). The top panel contains the realistic case 
with proper admixtures of both color-singlet and -octet contributions, and shows a strong J/ij} component, yet being 
dominated by Xc productions with a ratio P{J/'4’)/PiXc) ~ 2. The middle (bottom) panel exhibits the contributions 
from initial QQ color-singlet (color-octet) pairs alone. The patterns are strikingly different, as is already visible in 
the earlier Fig. 4. Golor-singlet production of the initial QQ results to largest extent in Xc-formation (mostly Xci), 
whereas T/'0) Vc and he are down by far more than an order of magnitude. Gontrary to that, when the initial QQ is 
produced in a color-octet, the situation is almost reversed. J/^/> production is the largest among the individual states, 
and is about equal to the sum of all three Xc states. 

Figs. 10 and 11 detail the beam energy (-^/s) dependence of color-singlet and -octet contributions for the charmonium 
and bottonium sector, respectively. In Fig. 10 the relative probabilities for Xc and J/'0 production are shown in the 
top panel, whereas the actual cross-sections are plotted in the bottom panel — both against the pp {pp) center-of-mass 
energy Ecm = Vs. Gorrespondingly, Fig. 11 depicts the behavior of Xb versus T production. The fat lines are the 
sum of color-singlet and color-octet contributions, the thin dashed (dotted) lines represent the color-singlet (-octet) 
contributions alone. From Fig. 10 one reads off that, in the present model, the color-singlet component of J/'0 is 
almost constant with if cm, whereas the color-octet contribution rises strongly with increasing if cm- At the same time 
color-singlet and -octet contributions to Xc are of the same order of magnitude, and decrease with ifem- As a result, 
at ISR energies 30 - 63 GeV, the J/V’ production is highly suppressed as compared to the Xc yield, at Tevatron energy 
1.8 TeV, the two channels become almost comparable, and at LHG energy 14 TeV, they become about equal. A 
similar pattern is seen in Fig. 11 for Xb and T, however with an even more pronounced energy dependence. Again, Xb 
production is by far favored over T formation at lower energies, but now already at Tevatron energy the two states 
are equally populated, and at LHG energy the T production even dominates by a factor of ~ 1.5. 
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FIG. 9. Relative proportions of the various included charmonium states at Tevatron energy = 1.8 TeV: top panel 
contains the realistic case with jrroper admixtures of both color-singlet and color-octet contributions, middle panel exhibite 
the contributions from initial QQ color-singlet pairs alone, whereas bottom panel displays the contributions from initial QQ 
color-octet pairs only. 
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FIG. 10. Relative probabilities for Xc and J/'il; production (top), and corresponding cross-sections (bottom) versus pp (pp) 
center-of-mass energy Ecm = \/s- The fat lines are the sum of color-singlet and color-octet contributions, the thin dashed 
(dotted) lines represent the color-singlet (color-octet) contributions alone. 
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FIG. 11. Relative probabilities for Xb s-ncl T production (top), and corresponding cross-sections (bottom) versus pp (pp) 
center-of-mass energy Ecm = sfs- The fat lines are the sum of color-singlet and color-octet contributions, the thin dashed 
(dotted) lines represent the color-singlet (color-octet) contributions alone. 
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C. Importance of parton evolution and gluon radiation 


Next, I turn to the second question raised above which concerns the importance and effects of gluon radiation 
during the PQCD evolution. As illustrated in Fig.2, gluon bremsstrahlung from both the virtual gluon (in case of the 
gluon fragmentation process) and from the produced QQ pair are included within the standard DGLAP evolution plus 
angular-ordering in the leading log approximation ||^ . This (multi-)gluon emission is a most important aspect in the 
present approach. Firstly, it changes the color-state of the QQ which makes it possible for an initial QQ color-octet 
to evolve into a final QQ color-singlet and to coalesce to a physical bound-state. Without gluon emission, only 
those QQ pairs that are initially produced as color-singlets could contribute to quarkonium formation. Secondly, 
gluon radiation causes significant energy loss of the QQ system at high energy, causing the relative momentum of 
the pair to be smaller in the mean, and their relative distance to increase slower. In effect, this dynamical feature 
enhances the probability for coalescence, since configurations with very small relative momentum and separation are 
strongly favored (c.f. Fig. 6). In other words, the presence of gluon radiation is crucial, because it allows the natural 
inclusion of both color-singlet and color-octet initial QQ pairs to evolve into quarkonium, and because the coalescence 
of final QQ systems is more likely to occur after sufficient gluon radiation to reduce the relative velocity to a value 
commensurate with the non-relativistic kinematics of the accessible quarkonium bound-states. 

Fig. 12 clearly displays the effect of gluon radiation in the model calculations. Confronted are the production 
cross-sections (^) for the realistic scenario including gluon emission during the PQCD evolution which are the same 
as in Fig. 7 (fat lines), with the corresponding cross-sections obtained by switching off all gluon radiation (thin lines). 
The comparison of the calculations with and without multiple gluon emission reveals two substantial differences. First, 
the overall magnitude in the case without gluon radiation is significantly lower by a factor of about 3 for J/'ij) and 
T, but only slightly less for Xc- Second, the shapes of the S-wave states J/'ij} and T are quite different in the two 
scenarios, indicating a softer, i.e. steeper, pj_-dependence in the no-radiation case. The shape of the cross-section for 
the P-wave Xc, on the other hand, shows a suppression at low pj_ in the case without gluon radiation, but is hardly 
affected at large py. 

Both these points can be intuitively understood. The substantial under-population of J/ij} and T in the no¬ 
radiation case is due to the fact that without gluon emision, only the color-singlet component of initial QQ production 
contributes, because color-octet pairs have no chance to convert into a final-state color-singlet, as is necessary to form 
quarkonium. Since at ^/s = 1.8 TeV the color-singlet cross-sections are smaller than the color-octet counterparts (c.f 
Figs. 10 and 11), the yield of J/ij} and T is naturally underestimated. On the other hand, the Xc production is much 
less deviating in the two cases, because here color-singlet and color-octet cross-sections are roughly of equal magnitude. 

The main lesson learned from Fig. 12 is that the neglect of gluon radiation results in quarkonium cross-sections 
which are incompatible with the data shown earlier in Fig. 7, both with respect to the absolute yield and the shape 
of the py-dependence. 

Fig. 13 further confirms the importance of multi-gluon emission. Plotted are the probability distributions of the 
number of gluons being emitted during the PQCD evolution between the hard vertex ab —>■ c*d of Fig. 2, up to the 
point of coalescence to a quarkonium bound-state. The full (dashed) histograms are for charmonium (bottonium) at 
pp (pp center-of-mass energies 63 CeV (top), 1.8 TeV (middle), and 14 TeV (bottom). One observes that bb evolution 
is generally accompanied by a larger amount of gluon radiation than cc, an effect that increases with energy. Even at 
ISR energy 63 CeV, gluon radiation is non-negligible, with mean values of (Ng) = 1.4 (1.7) for cc {bb). At Tevatron 
energy 1.8 TeV, the mean number of gluons emitted per QQ system are (Ng) = 1.9 (4.7) for cc {bb), and at LHC energy 
14 TeV, one has {Ng) = 2.1 (5.8). This shows that even at ISR energy 63 CeV, gluon radiation is non-negligible, and 
that it is most essential at Tevatron and LHC. 


** Notice that a color-octet QQ produced in a state must radiate at least two gluons in order to convert into a color-singlet 
state with quantum numbers ot Jftp or T. On the other hand, if it is produced in a 3j state, the emission of only one gluon is 
necssary to form a XcJ or Xbj- 
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FIG. 12. Quarkonium production cross-sections of the previous Fig. 7 revisited: Fat lines represent the realistic scenario 
including gluon emission during the PQCD evolution, whereas thin lines correspond to the no-radiation scenario with all gluon 
radiation switched off. Top part includes prompt J/ip’s only, middle part shows Xc’s from radiative Xc ^ 7 + and bottom 
part contains T’s, including the radiative feed-down from the Xb states. 
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FIG. 13. Probability distributions of the average number of gluons being emitted during the PQCD evolution between the 
hard vertex ab c*d of Fig. 2, up to the point of coalescence to a quarkonium bound-state. The full (dashed) histograms are 
for charmonium (bottonium) at pp {pp center-of-mass energies 63 GeV (top), 1.8 TeV (middle), and 14 TeV (bottom). 
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D. Characteristics of the space-time development 


I finally address the third and last question raised at the end of Sec. III.A, namely, how does the time evolution of 
the QQ system proceed after being initially produced either directly at the hard vertex (Fig. a), or by virtual gluon 
decay (Fig. 2b), as it propagates in space up to the point where it either forms a bound-state, or ends up in other 
hadrons as open charm, respectively bottom. In view of the preceding discussions, I wish to detail some characteristic 
space-time features that are related to the issues of multi-gluon emission and color-singlet/color-octet processes. 


E = 1.8 TeV 

cm 




FIG. 14. Time-evolution of the average relative QQ-momentum = (|kQ — k^|) (top) and of the average QQ-separation 
Cgg = (|rQ — r^l) (bottom). Solid curves represent the ‘full QCD evolution’ including gluon radiation, and dashed curves 
correspond to_‘free propagation’ without gluon radiation. The plots refer to charmonium production at 1.8 TeV and include 
only those QQ that actually end up as hnal charmionium states. 
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Fig. 14 vividly illustrates the different time-evolution patterns for the ‘full QCD evolution’ including gluon radia¬ 
tion (solid curves) and the ‘free propagation’ without gluon radiation (dashed curves), here exemplarily depicted for 
charmonium production at Tevatron energy 1.8 TeV. Only those QQ are included in the plots that actually end up 
as final charmonium states, i.e., open charm production has been subtracted. The top panel shows the change of 
the average relative momentum fcgg = (|kg — k^|) of Q and Q, as a function of time in the hadronic center-of-mass 

frame. The bottom panel displays the corresponding time-dependence of the relative separation of Q and Q in the 
pair rest-frame, rgg = (|rg — rg-|). The two plots confirm what was mentioned already ealier, namely that multiple 

gluon emission leads to a significant energy loss which results in a drop of the relative momentum fcgg (transverse 
momentum broadening is comparably ineffective). The natural consequence is a slower time change in the relative 
pair distance rgg in the case with gluon radiation. Without the emission of gluons, the Q and Q do not alter their 
initially obtained momentum and propagate classically in free space, which leads to a faster growth in their spatial 
relative separation. 

This behavior has interesting consequences for the coalescence process and the resulting quarkonium yield. As is 
evident from the previous Fig. 6, QQ configurations with small rgg and small /cgg have the largest coalescence 
probabilities, and hence are strongly preferred in the formation of quarkonium bound-states. This in turn explains 
why in Fig. 12 the no-radiation scenario gives a lower quarkonium cross-section than the case with gluon radiation. 

Fig. 15 exhibits, again for the charmonium sector at 1.8 TeV, further striking differences between the ‘full QCD 
evolution’ including gluon radiation, and the ‘free propagation’ without gluon emission. The top part shows the time 
rate of change of the probability for QQ bound-state formation, versus the probability for Q and Q to disappear 
individually in open charm (i.e., in other charmed hadrons by recombination with partons from the rest of the pp 
event). The bottom part depicts for those QQ pairs that indeed form bound-states, the time-dependence of the 
individual coalescence probabilities (n^®)(t)) for the spectrum of charmonium states i Sq,^ S'l,^ Pi,^ Pj. Both plots 
exhibit very different shapes for the two cases, with and without gluon radiation. Although this time-development is 
not directly accessible in experiments, one may think of sensible observables that reflect some of these features and 
would allow to study the space-time structure by indirect measurement. 

Finally, Fig. 16 gives an instructive summary of the space-time history that the initially produced QQ pairs have 
experienced on their way to coalescence as final QQ color-singlet systems. Again, I choose as example the charmonium 
sector at pp center-of-mass energy of 1.8 TeV. The left side shows the distributions in relative separation ri„i, relative 
momentum kini and invariant mass Mini of the QQ pair upon production, and the right side displays the corresponding 
spectra with respect to rfin, kfm and Mfin of the final QQ upon coalescence to charmonium bound-states. 

The plots are rather self-evident and in line with intuitive expectations: The initial Q and Q are highly localized 
with a mean spatial extent of ~ 10“2 /m, whereas the final QQ pair has a much larger typical size, being of the order 
of the size of the charmonium states (c.f. Table 1). On the other hand, the initial and final spectra in the relative 
QQ momentum show a reverse pattern, namely a significantly narrower distibution of the final QQ momentum upon 
coalescence, being the cumulate effect of gluon radiation that decreases that decreases the originally obtained relative 
momentum of the QQ when it was produced initially. Similarly, the mass spectrum of the initial QQ, being related 
to the distribution of Q^ values upon production, spreads out to masses well above 10 GeV, although exponentially 
damped. The mass distribution of the final QQ, on the other hand, is concentrated in its mean around the 3-4 GeV 
region, where the charmonium states lie. 


It should be noted that the pronounced curvature of around 1 fm arises as a consequence of the formation time of 
emitted gluons, tform oc Eg/k\g. Since the bulk of bremsstrahlung is associated with low-energy gluons, the cumulate effect of 
these latter on fegg becomes significant only after the typical formation time 1 fm. 
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FIG. 15. Characteristics of the space-time evolution of the ‘full QCD evolution’ including gluon radiation, and the ‘free 
propagation’ without gluon emission for charm production at pp center-of-mass energy o_^1.8 TeV. Top part: time-development 
of the probability for QQ bound-state formation, versus the probability for Q and Q to hadronize in open charm. Bot¬ 
tom part: time-dependence of the individual coalescence probabilities {11^*^ (t)) for the spectrum of charmonium states 
i So,^ Si,^ Pi,® Pj. Included are only those QQ pairs that indeed form bound-states, open charm production is subtracted. 
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FIG. 16. Characteristics of the space-time evolution reflected in initial and final QQ spectra, for the case of charmonium 
pp center-of-mass energy of 1.8 TeV. Left panel: Distributions in relative separation Vini = |rQ — relative momentum 

kini = I kg — k^|t=to and invariant mass Mini = {kQ of the initial QQ pair upon production at t = to. Right 

panel: Distributions in the corresponding quantities Vfin, kfin, Mfin of the final QQ pair upon coalescence. 
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IV. CONCLUDING REMARKS AN OUTLOOK 


A. Summary 


The scope of this work_concerned the development and application of a space-time description for the production of 
heavy quarkonium as QQ bound-states in hadron-hadron collisions. The theoretical framework was developed on the 
basis of the common factorization assumption of short-distance PQCD production and evolution of heavy QQ pairs 
and the long-distance bound-state formation process leading to final-state quarkonia. 

The separation of perturbative and non-perturbative momentum scales (and space-time scales) was exploited to 
formulate the dynamics of quarkonium production and the associated cross-section as the product of three functions 
corresponding to QQ production, parton cascade evolution, and coalescence, repectively. The practical calulation 
involved three corresponding elements: hrstly, the ‘hard’ production of the initial QQ, directly or via gluon fragmenta¬ 
tion and including both color-singlet and -octet contributions, secondly, the subsequent development of the QQ system 
within a space-time generalization of the DGLAP parton-evolution formalism in position- and momentum-space, and 
thirdly, the actual formation of the bound-states being described by the space-time dependent overlap between the 
QQ pair and the spectrum of quarkonium wave-functions. 

This approach was then applied to pp (pp) collisions in the energy range ^/s = 30 GeV - 14 TeV, corresponding 
to the range of center-of-mass energies between ISR, Tevatron, and LHC experiments, by simulating the space-time 
development of the complete particle system of a collision event, using a suitably extended version of the computer 
program VNI. Comparison with experimental data from ISR and Tevatron for both charmonium and bottonium 
production revealed very good agreement of the model calculations, essentially without the necessity of introducing 
tunable, free parameters. Details of the underlying dynamics that lead to the decent postdictions of ISR and Tevatron 
data (and predictions for LHC), were analyzed in particular with respect to the relevance of gluon radiation, the 
contributions of color-singlet and color-octet QQ systems, and the characteristics of the space-time structure. 


B. Future perspectives 


Finally some remarks on future perspectives. The good agreement with available experimental data of the present 
approach to quarkonium production in pp or pp collisions, as well as its transparent physics concept, suggests to extent 
applications and tests of the model to other collision systems. 

• pion-proton collisions are experimentally well studied too in terms of quarkonium production (although at lower 
energies), observed with even distinctly larger yields, partly due to the different valence quark content of the 
initial state. It would be interesting to test whether the model advertized here is capable of describing these 
reactions equally well. 


• proton-nucleus collisions are another extension for which again quite some experimental data exist for various 
energies and different nuclei. The novel feature here, as compared to hadron-hadron collisions, is the propagation 
of the QQ pair in nuclear matter, that is, before forming a bound-state, Q and/or Q may experience scatterings off 
partons belonging to close-by nucleons in the traversed nucleus. As a consequence, one expects the quarkonium 
yield to be reduced due to the deflection of Q and/or Q and the resulting failure of coalescence. 


• nucleus-nucleus collisions, finally, are a most interesting ground for studying quarkonium production in QCD 
matter, in particular since it has become appreciated as one of the best experimental probes for indirect detection 
of a quark-gluon plasma, a high-density and high-temperature state of deconfined partonic matter for which first 
hopeful signals have been extracted recently in lead-on-lead collisions at the CERN SPS facility. Compared to 
proton-nucleus collisions, in heavy-ion reactions several new physics issues arise. For example, the large size of the 
nuclear collision system makes it likely that the traversed distance of a heavy QQ system can extend over several 
fm so that not only the Q and Q are likely to re-interact before they may coalesce, but also the coalesced QQ 
bound-states can re-collide and being broken up again (so-called ‘dissociation’). Moreover, in central collisions 
of heavy nuclei at sufficiently large energy per nucleon, it is expected 241 that copious parton materialization 
and mini-jet production through very frequent parton collisions, plus further parton multiplication through 
gluon emission induced by partonic re-scatterings, may create very quickly (with in a fraction of a fm) a very 
dense and highly excited partonic environment. Whether or not this QCD matter approaches an equilibrium 
state, i.e. some sort of quark-gluon plasma, the properties of the hot and dense parton matter are of their 
own fundamental interest. In this context, quarkonium production is certainly a pre-destined probe to study 
these properties, as mentioned before. Both the yield of formed QQ bound-states and the survival probability 
of resulting quarkonium states are likely to be suppressed, the degree of suppression being an indicator of the 
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density and temperature of the nuclear matter. It is clear that here the space-time picture is an inevitable 
necessity to study quarkonium production in high-density QCD, where accounting for the dynamically evolving 
space-time structure of the collisions is most important, and where the time-dependent spatial profile of the 
particle density is the key quantity that determines the fate of quarkonium. Since the local density is created 
by the internetted space-time structure of many parton cascades evolving close-by in the central collision region, 
it the only way of ‘calculating’ the dynamics is by rigorous accounting for both the space-time variables and 
the energy-momentum variables for each particle present at a given time during the collision. Here lies a real 
challenge for the model advocated in this paper. The technicalities of extending the present status to nuclear 
collisions are straightforward, though non-trivial, and work is in progress. 
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